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Abstract 

An attempt is made to generalise the ideas introduced by Haldane and others regarding Bosonizing the 
Fermi surface. The present attempt involves introduction of Bose fields that correspond to displacements of 
the Fermi sea rather than just the Fermi surface. This enables the study of short wavelength fluctuations of the 
Fermi surface and hence the dispersion of single particle excitations with high energy. The number conserving 
product of two Fermi fields is represented as a simple combination of these Bose fields. It is shown that most(!) 
commutation rules involving these number conserving products are reproduced exactly, as are the dynamical 
correlation functions of the free theory. Also the work of Sharp, Menikoff and Goldin has shown that the field 
operator may be viewed as a Unitary representation of the local current group. An explicit realisation of this 
unitary representation is given in terms of canonical conjugate of the density operator. 

1 Introduction 

Recent years have seen remarkable developments in many-body theory in the form an assortment of techniques 
that may be loosely termed bosonization. The beginnings of these types of techiniques may be traced back to 
the work of Tomonaga[[l] and later on by Luttinger Q and by Leib and Mattis ||. In the 70's an attempt was 
made by Luther || at generalising these ideas to higher dimensions. Closely related to this is work by Sharp et. 
al. in current algebra. Mention must also be made about the work of Feenberg and his collaborators |Q who 
have a theory of correlated electrons written down in the standard wavefunction approach familiar in elementary 
quantum mechanics. It seems that the approach adopted here has some similarites with this work although the 
details are not identical. More progress was made by Haldane || which culminated in the explicit computation 
of the single particle propagator by Castro-Neto and Fradkin || and by Houghton, Marston et.al. @ and also by 
Kopietz et. al. J8|. Rigorous work by Frohlich and Marchetti|J is also along similar lines. 

The attempt made here is to generalise the concepts of Haldane [|J to accomodate short wavelengths fluctu- 
ations where the concept of a linerised bare fermion energy dispersion is no longer valid. To motivate progress 
in this direction one must first introduce the concept of the canonical conjugate of the Fermi density distribu- 
tion. This concept is likely to be important in the construction of the field operator in terms of the Fermi sea 
displacements, although the latter is not completed in this article. The concept of the velocity operator being 
the canonical conjugate of the density has been around for a long time, and this has been exploited in the study 
of Hell by Sunakawa et. al. [ fTpf . However, the author is not aware of a rigorous study of the meaning of this 
object, in particular, an explicit formula for the canonical conjugate of the density operator has to the best of my 
knowledge never been written down in terms of the field operators. The work by Sharp et. al. JTT| comes close 
to what I am attempting here. 



2 Unitary Representation of the Local Current Group 

The work of Sharp et. al. [jlTJ established that the current algebra itself does not convey the underlying particle 
statistics rather the statistics is hidden in a choice of a unitary representation of the local current group. Here, I 
try to write down some formulas that provide (possibly for the first time) an explicit and exact representation of 
the Fermi field operator in terms of the canonical conjugate of the density operator thereby providing an explicit 
realisation of the claims of Sharp et. al. UM 
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2.1 Some Mathematical Identities 

is a bose-like object; namely it satisfies 



Notation: 



[p(xa),p(x a )] = 
y dx p(xcr) 



n a = — 

Pa = (n a ) 

Pa is a c-number, whereas N a and are operators. 

So it natural to introduce the canonical conjugate n(xcr). 

[n(xa),p(xV)]= l( 5 d (x-xX y 

[n(xo-),n(xV)] = o 

Propose the ansatz 

?/>(xer) = exp(— Hl(x.a))exp(iQ([p\: xct))(/9(xct)) = 
The above Eq.( ||) shall be called the DPVA ansatz. Here again, 

Il(xcr) = i ln[(p3 + (S^(xcr))(p CT + Sp(xa)) 5 exp(-i$([p]; xcr))] 

Sip(ynj) = ip(xa) - 
5p(xcr) = p(xcr) - 

p CT is a c-number given by p a = (n a ). Where $([/>]; xcr) is some hermitian functional to be computed later, 
ansatz automatically satsfies 

(XCT)V'(XO') = yO(xcr) 



and 
Write 



where n a = and 



[V>(x<j),p(xV)] = 5 d (x-x')5 a y^(xa) 
p(xo-) = n CT + — ^ Pqaexp(-iq.x) 



q#0 



Therefore 



where 



n(xcr) = X 0cr + exp(iq.x)X qcr 
[X 0CT ,iV CT ,]=z^ 

^(xcr) = esp(-«^eafp(«q.x)X qff )F([{/9 kcr /}];xcr) 



= e^(^([{p kff '}];xa))(n CT + 1 ^ 



exp^-iq.xjpq^j 



q#0 



It is possible to write these formulas in momentum space. This will be useful later on. First we write 

*([p]; xcr ) = ^Hipli^^pi-i^) 



Therefore 

tp(x<r) = exp(-i^exp(ici.x)Xq a )exp(i^exp(-ici.x)(j)([p}-qa))(n a + — ^ pq a exp(-iq.x))? (18) 
q q q/o 

It is possible to write the corresponding formula in momentum space by making the identifications 

exp^.x) - T_ q (k) (19) 

exp(-iq.x) -» T q (k) (20) 

where 

T q (k) - ea:p(q.V k ) (21) 

V-(ka) = e^(-^T_ q (k)X qCT ) ea ;p(^T q (k)^([p];q ( 7))(Ar (7 + ^p qCT T q (k))^ M (22) 
q q q#o 

The translation operators translate the Kronecker delta that appears in the extreme right. It may be verified 
that 

[V(k<7),p q ^]=V(k- q( 7)^ (23) 
M)«T),N a ,}=iP(k*)S aia , (24) 

Here 

^(xcr) = — r V ea;p(ik.x)V'(kCT) (25) 
^ k 

Also the following identities are going to be important. 

^(ka) = ezpH^T_ q (k)AV)./W(M) (26) 
q 

f ka ([p\) = exp(zJ2W<P([p};qcT))(N a + ]T T q (k)p qff )i<J k , (27) 
q q#o 

Check by expansion that 

E/k +q /2.(M)./k- q /2.(M) - Pq. (28) 
k 

2.2 Proof (??) of the Fermion Commutation Rules 

We use the theory of distributions to ( try ) prove rigorously (good enough for a physicist!) the fermion anticom- 
mutation rules. Pardon the pretense of mathematical rigor. The approach is as follows. Let be the space of all 
smooth functions from S — R 3 x {T,|} to C. Further, on this space define the inner product (Schwartz space). 
In : L 2 C (S) x L 2 C (S) -» C 



(f\g) 

Since the Fermi fields are operator-valued distributions we can construct operators 

C f = / dx X^(xcr)/(xcr) 

and we then have to prove 

{c f ,c g } = (29) 
c) = (30) 
The above two relations shall be called c — c anticommutation rules. 

{c /)C t} = <<?!/> (31) 
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The above relation shall be called the c — anticommutation rule, for all /, g that belong to Lq(S). The claim 
is, that these relations are satisfied provided $ obeys the recursion below. 

*([{p(yi^)-«(yi-x')* ai y}];»7) 

+$([/9];xV) - $([/o];xct) 

- ^([{p(yio-i) - *(yi -*)<W}];xV) ^tott (32) 

where m is an odd integer. This recursion is to be satisfied for all (xcr) ^ (x cr ). At precisely (xcr) = (x cr ) 
the recursion obviously breaks down. But this is not a serious drawback as we shall soon find out. The proof 
involves working with a basis. Here the integer m has to be even for bosons and odd for fermions. The imposition 
of the commutation rules by themselves do not provide us with a formula for the hermitian functional Any 
redefinition of II (xcr) consistent with [II, p] = iS(...) may be absorbed by a suitable redefinition of <!> (akin to 
gauge transformations pointed out by AHC Neto, private communication). Solution to $([p];qcr) by making 
contact with the free theory. A random choice of $([p];xcr) that satisfies the recursion is (A. J. Leggett : private 
communication) the Leggett ansatz (for fermions) (Also F.D.M. Haldane 

cr 

redefinitions of $ that leave the statistics invariant are the 'generalised guage transformations'. Disclaimer:Please 
pardon my pretensions at rigor in the following paragraph. It is easy to get carried away. The relation that 
relates two fermi fields by a generalised guage transformation is an equivalence relation. Therefore , analogous to 
the claim that a state of a system is a ray in Hilbert space, one is tempted to call the equivalence class of fermi 
fields under these transformations the fermi distribution. Therefore a fermi distribution is not just one fermi 
field but a whole bunch of equivalent ones. Therefore a state containing one fermion at a space point is obtained 
by acting the creation fermi distribution on the vacuum. Thus one may construct the Fock space by repeatedly 
acting these fermi distributions on the vacuum. Therefore there are two ways by which one can choose $. One 
is a random choice which satisfies the recursion. In which case LT is determined by Eq.(^|). The other choice is 
II = In which case, $ can no longer be chosen arbitrarily. It has to be determined by making contact with 
the free theory. 

Here 9{x) is the Heaviside step function. The unpleasentness caused by the fact that the recursion for fermions 
does not hold when (xcr) = (x a ) is probably remediable by multiplying the fermi fields which are operator- valued 
distributions by arbitrary smooth functions and proving properties about these latter objects. More importantly, t 
has to invertible in order for the ansatz to satisfy the recursion. That such a bijective mapping exists is guaranteed 
by the theory of cardinals. However, this mapping is not continuous let alone differentiable and therefore of little 
practical value. (A continuous mapping would imply a homeomorphism between R and R 3 eg.) 

The claim is that the state of the fermi system is prescribed by prescribing an amplitude for finding the system 
in a given configuration of densities. Let 

Wfs{\p]) = ®H({p])exp(-U F s[p])exp(i6 F s([p])) 

be the wavefunctional of the noninteracting fermi sea. The Ufs[p] is uniquely determined by prescribing all 
moments of the density operator. The functional 0fs[p] cannot be so determined. Moreover, since the amplitude 
for finding the system with negative densities is zero, we must also have a prefactor, 

e H ([p}) = n yiai o H {p(yi<ri)) 

Where 9h(x) is the Heaviside unit step function. 

2.3 Making Contact With the Free Theory 

Consider the operator 

n(ya) = / dx ^ f (x + y/2a)ip(x - y/2a) = ip^ (ka)^(ka)exp(-ik.y) (33) 

J k 
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n(ya)W FS ([p\) = n (ya)W FS ([p}) (34) 
^o(yo-) = ^2 0{k F - k)exp{-i\t.y) (35) 

k 

First assume y ^ 0. In eq. ( |34| ) there are two undetermined functionals. One is $([/>]; xer) and the other is 
0fs([p])- We have remarked earlier that redefinitions of <&([/>]; xer) are similar to 'gauge transformations'. These 
functional gauge transformations leave the local density of particles invariant but alter the statistics. There is a 
subgroup of these functional gauge transformations that also leaves the statistics invariant. A realisation of this 
subgroup is achieved by for example: 

mp); X a)^m P };xa)+e([p}) 

This may be exploited to our advantage while solving eq.( |34|). In other words 9 F s([p\) can be absorbed into 
$([p]; xer) without altering the statistics. This leaves us with an equation for $([/>]; xer) in terms of the properties 
of the free theory exclusively. Exploiting the recursion relation for 3>([p]; xer) we arrive at the result (Here y ^ 0.) 

J c?x (p((x + y/2V))4(p((x - y/2)a)f* exp[i${[p]; (x - y/2)a) - i$(\p]; (x + y/2)a)} 

exp(U FS [{p(y'a) - S(y - (x - y/2))^}] - U FS [{ P (ya) - 5(y - (x + y/2))^, JD 
e([{p(y'a)-5(y -( X + y/2))S a/ J}) 

WRT^MJ = ~ Mya) (36) 

The terms involving the Heaviside are indeterminate unless we expand around 

p(ya) = 6 d (0) + Sp(ya) (37) 

and, 

( p(ya) ) = S d (0) = p u (38) 



2.4 Point Splitting or no Point Splitting ? 

The attempts made here are partly based on the work of Ligouri and Mintchev on Generalised statistics |H| and 
work of Goldin et. al.|n]] and the series by Reed and Simon[jl5| Here I shall attempt to provide a framework 
within which questions such as the existence of the canonical conjugate of the fermi-density distribution may be 
addressed. For reasons of clarity we shall not insist on utmost generality. The philosophy being that the quest 
for utmost genarality should not cloud the underlying basic principles. We start off with some preliminaries. Let 
TL be an infinite dimensional separable Hilbert Space. We know from textbooks that such a space possesses a 
countable orthonormal basis B — {wi\i £ Z} . Here, Z is the set of all integers. ThusH = Set of all linear 
combinations of vectors chosen from B. We construct the tensor product of two such spaces 

This is defined to be the dual space of the space of all bilinear forms on the direct sum. In plain English this 
means something like this. Let / S 7i and g £ Tt define the object f®g to be that object which acts as shown 
below. Let < v, w > be an element of the Cartesian product TC x 7i. 

f(£)g < v, w >= (/, v)(g, w) 
Here, (/, v) stands for the inner product of / and v. Define also the inner product of two f(£)g and / §Qg 

= (/,/')(<?,</) 

Construct the space of all finite linear combinations of objects such as f®g with different choices for / and g. 
Lump them all into a set. You get a vector space. It is still not the vector space Because the space of 

all finite linear combinations of objects such as f®g is not complete. Not every Cauchy sequence converges. 
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Complete the space by appending the limit points of all Cauchy sequences from the space of all finite linear 
combinations of vectors of the type f®g. This complete space is the Hilbcrt space Similarly one can 

construct for n = 0, 1, 2, 3, ... Where we have set Tt° — C the set of complex numbers. Define the Fock 

Space over Ji as 

Physically, each element of it is an ordered collection of wavefunctions with different number of particles 

((po, <pi(ari), Lp 2 (xi,x 2 ), ip n , •••) 

is a typical element of T(TL). This is the Hilbert Space which we shall be working with. Let V n be the space of 
all decomposable vectors. 

P n = {/i(g)-(g)/n;/ieW} 



For each / eH define 



defined by 



b(f) :V n ^V n -\n> 1 
&*(/) :D" ->V n+1 ,n>0 

b(f) /i0...(g)/„ = V^(/,/i)/ 2 (g)...(g)/ n 



b*(f) /l(g)..-0/n = V^+l/(g)/l(g)/2(g)...(g)/„ 

We also define b(f)7i° = 0. By linearity we can extend the definitions to the space of all finite linear combinations 
of elements of V n namely C{V n ). For any ip e £(£>") and ^ e C{V n+1 ) 

II K/VII<v^ll/lllkll 



II &*(/)¥> ||< v^Tl || / mi v» || 
(^6*(/V) = (K/)V,v) 

So long as || / ||< oo, and &*(/) are bounded operators. An operator O is said to be bounded if 

sup M=1 || Oip\\ < oo 
is unbounded otherwise. The norm of a bounded operator is defined as 

II o ||= sup M =i || evil 

In order to describe fermions, it is necessary to construct orthogonal projectors on T(TL). In what follows c(/) 
will denote a fermi annhilation operator. c*(/) will denote a fermi creation operator. Physically, and naively 
speaking, these are the fermi operators in "momentum space" Ck and c£. First define P_ to be the projection 
operator that projects out only the antisymmetric parts of many body wavefunctions. For example, 

P-/i0/ 2 = i(A(g)/ 2 - / 2 (g)/l) 

We now have 

c(f) = P-b(f)P- 

C*(f)=P-b*(f)P- 

Let us take a more complicated example. Let us find out how c*(f)c(g) acts on a vector v = /i®/2- 

c*(f)c(g) = P_b*(f)P_P_b(g)P_ 

c*{f)c{g) = P.b*{f)P.b{g)P. 
c*(f)c(g)v = P_b*(f)P-b(g)P-v 
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^=^j(/l(8)/2-/2(g)/l) 

b{g)P_v = ^V2((g, h)f 2 - (g, f 2 )h) 

P_b(g)P_v = ±s/2{{g,h)f 2 - (g,f 2 )fi) 
b*{f)P.b{g)P.v = ^V2 2 ((g,/i)/(g)/2 - (<7,/ 2 )/(g)/i) 

c*(f) c(g)v = (^) 2 V2 2 ((.9, /i)[/(g)/ 2 - / 2 (g)/] - (.9, M/(g)/i - /i(g)/]) 

Having had a feel for how the fermi operators behave, we are now equipped to pose some more pertinent questions. 
Choose a basis 

B = {wi\i e 2} 

2.5 Definition of the Fermi Density Distribution 

Here we would like to capture the notion of the fermi density operator. Physicists call it p(x) — ip*(x)ip(x). 
Multiplication of two fermi fields at the same point is a tricky business and we would like to make more sense out 
of it. For this we have to set our single particle Hilbert Space: 

U = L p (TZ 3 )0W 

Here, L P (1Z 3 ) is the space of all periodic functions with period L in each space direction. That is if u e L P (1Z 3 ) 
then 

u(xi + L,x 2 ,x 3 ) = u(xi,x 2 ,x 3 ) 
u(xi,x 2 + L,x 3 ) = u(xi,x 2 ,x 3 ) 
u(x 1 ,x 2 ,x 3 + L) = u(x 1 ,x 2 ,x 3 ) 
W is the spin space spanned by two vectors. An orthonormal basis for W 

A typical element of H is given by /(x) ® A basis for Tt is given by 

T 



B = {\J ^3 ea; p(*qn-x)(g)6;n = (ni,n 2 ,n 3 ) e Z 3 ,s e {?,!}; 

,27ml 2nn 2 2nn 3 

We move on to the definition of the fermi-density distribution. The Hilbert Space is the space of all 

n-particle wavefunctions with no symmetry restrictions. From this we may construct orthogonal subspaces 

H® n = P + H® n 

Tensors from are orthogonal to tensors from H_ . The only exceptions are when n = or n = 1. 



nf Q =nf Q =c 
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The space W®"is the space of bosonic-wavefunctions and the space H® n is the space of fermionic wavefunctions. 
The definition of the fermi density distribution proceeds as follows. Let v be written as 

v = a ( cr )^ 

<re{T,i} 

The Fermi density distibution is an operator on the Fock Space, given a vector / v G H in the single particle 
Hilbcrt Space, and a tensor tp in the n-particle subspace of ofF(H), there exists a corresponding operator p(f v) 
that acts as follows: 

[p(.f v)(p\ n (x 1 <j 1 ,x 2 V2, ....,x n cr„) = 

if ipe U® n and 

n 

, X2*T2, , X n (T„) = ^ / (Xj ) q (<Tj ) y n (Xj C7i , x 2 Q-2 , ....,X n <7„) 

1=1 

when e . The physical meaning of this abstract operator will become clear in the next subsection. 

2.6 Definition of the Canonical Conjugate of the Fermi Density 

We introduce some notation. Let g = y^-exp(ik m .x) ® £ r 

ip(k m r) = c(g) 
p(k m r) = p(g) 

This p(k m r) is nothing but the density operator in momentum space, familiar to Physicists 

We want to define the canonical conjugate of the density operator as an operator that maps the Fock space(or a 
subset thereof) on to itself. If <p e Tt° = C then 

= 

Let ip e n = 2,3, ... then 

Xq m ^ = 

The important cases are when p e H™, n = 2, 3, ... or if ip e W. In such a case, we set n = N® ^ for s e {J., |}. 
Let us introduce some more notation. N s = Pq = o s is the number operator to be distinguished from the c-number 
N®. The eigenvalue of N s is N® when it acts on a state such as p> € HI. Some more notation. 



5 V(k m s) = V(k m s) - y/NSSu^fi 

and 

(5 p(k m s) = p(k m s) - 7V s °<5 km ,o 
The Canonical Conjugate of the density distribution in real space denoted by n s (x) is defined as follows. 

II s (a;i + L,x 2 ,x 3 ) = U s (x 1 ,x 2 ,x 3 ) 

IT s (a;i,a;2 + L,x 3 ) = U s (x 1 ,x 2 ,x 3 ) 
U s (x 1 ,x 2 ,x 3 + L) = U s (x 1 ,x 2 ,x 3 ) 

The definition is as follows. 



8 



X qmS =Hn[(l + — ^=^5V(k„s)T_ kn (q m ))(l+^ 
^o5p(k nS )T kn (q m ))-3e a :p(-^0([ /9 ];k nS )T k „(q m ))]5 qm ,o (39) 



k„ 

$([p];xs) = y^^([jo];k n s)ea;p(-ik n .x) 

k„ 

7k n (qm) = ea;p(k n .V qm ) 

The translation operator translates the q m in the Kronecker delta that appears in the extreme right by k n and 
$([,o];xs) satisfies a recursion explained in detail in the previous manuscript. The logarithm is to be interpreted 
as an expansion around the leading term which is either N® or J N®. The question of existence of Xq mS now 
reduces to demonstrating that this operator (possibly unbounded) maps its domain of defmition(densely defined 
in Fock space) on to the Fock space. Defining the limit of the series expansion is likely to be the major bottleneck 
in demonstrating the existence of Xq mS . That this is the canonical conjugate of the density operator is not at 
all obvious from the above definition. A rigorous proof of that is also likely to be difficult. Considering that we 
arrived at this formula by first postulating the existence of II s (x), it is probably safe to just say " it is clear that 
" II s (x), is in fact the canonical conjugate of p. The way in which the above formula can be deduced may be 
motivated as follows: 

^(xc) = 77T 52 exp(ik.x)ip(ka) 

k 

= exp(—i 52 exp(iq.x.)Xq a )exp(i 52 ex P{~ i ( l-' x -)4>([p]'-, ( ¥ 7 )) 
q q 

( n<7 + v 52 Pqaexp(-iq.x))^ (40) 

In the above Eq.( |i(]) ONLY on the right side make the replacements, 

exp{iq.x) -> T_ q (k) (41) 
exp(-^q.x) -» T q (k) (42) 

where 

T q (k) = exp{q.V k ) (43) 

and append a 5 k o on the extreme right. Also on the LEFT side of Eqn.( ^) make the replacement ?/>(xcr) by 
i/)(ka). And you get a formula for ^(ker). In order to get a formula for Xq mS we have to invert the relation and 
obtain, 

n(xs) = i ^(v^vj + 52 

exp(ik n .x.)Sip(k n s)) 

k„ 

ea;p(-z52 ea:: P(-*k n .x)0([p];k n s))(7V° + 52 ^k„ s exp(-ik n .x))"5] (44) 

k„ k„ 

Make the replacements on the right side of Eqn.( |44|) 

exp(ik n .x) -> T_ kn (q m ) (45) 

expHk n .x)^T kn (q m ) (46) 

where and append a i5 qm! o on the extreme right. Also on the LEFT side of Eqn.( |44|) replace II(xs) by Xq mS . 
This results in formula given in Eqn. ( p9|). 



9 



3 Expression in terms of Fermi Sea Displacemnts 

The claim is that the following exact relation holds, 



where, 

and 

and, 



t N N + 
=k+q/2 c k- q /2 = rip (k) <5 q>0 + (W ^y)[A k (q)0(-k.q)a k (-q) + A k (-q)0(k.q)a k (q)] 

+ E E r^;^(k,q)a ki ( qi )a k2 (q 2 ) (47) 
ki,k 2 qi,q2 

{c k ,c k '}-0 (48) 

{c k ,c k ,} = <5 kik , (49) 

[a k (q),a k /(q')] =0 (50) 

K(q),a k '(q')] = 5 k,k'<W ( 51 ) 



and, 
and 



A k (q) = y/n F (\a + q/2)(l - n F (k - q/2)) (52) 
N = J2^W (53) 

k 

(A0=5> F (k) (54) 

k 

and n F (k) = #(fc/ — |k|) The 9(— k.q) is a reminder that A k (q) is nonzero only when k.q < 0. Also, 

[iV,a k (q)]=0 (55) 
[c p ,a k (q)] = something complicated (56) 

and the kinetic energy operator is, 

K = ^2 £kC k c k 

k 

in terms of the Bose fields, it is postulated to be, 

K = E o + E ^k(q)4(q)«k(q) (57) 

k,q 

w k (q) = A k (-q) — 
m 

and Eq — X) k £ k rl -F(k)- Also, the filled Fermi sea is identified with the Bose vacuum. \FS >— |0 >. The fact 
that this is true will become clearer later. At this point it is sufficient to note that w k (q) > for all k and q. We 
have to make sure that for the right choice of T, 

(Pq 1 (*l)Pq a (*2)P-q 1 -q a (*3)> + (58) 

The terms linear in the Fermi sea displacements in Eq.( f47j ) are chosen in order to ensure that the correct 
dynamical four-point functions are recovered. Now for the quadratic terms. The best way to derive formulae for 
them is to compute the exact dynamical six point function and set it equal to the expression got from the free 
theory. 

( c !+ q /2 W c k- q /2 (*)c k , +q , /2 (t )<v _«,' 11 (*' ) c l>> +q » /2 (*"K» _ q " /2 (*")) = J ( k > q. *; k ' > q' » *' ; k " . q" . *" ) 
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From the free theory we have 

"n / k.q x / k'.q' . k".q" » 

I (k, q, t; k , q , t : k ,q ,t ) — expii tjexpn 1 jexph 1 ) 

mm m 

[ (1 - n F (k - q/2))(l - n F {k - q'/2))n F (k + q/2)5 k+q/2!k »_ q » /2 5 k _q/2,k' + q'/2 ,5 k'-q'/2,k"+q"/2 

- ( 5k-q/2,k"+q"/2 5 k+q/2.k'-q72' 5 k'+q72,k"-q"/2( 1 - »f (k - q/2))n F (k + q'/2)nj?(k + q/2) ] (59) 

Also the six-point function may be evaluated in terms of the Bose fields(summation over superfluous indices is 
implied), 

J = (A k (q)a k (-q)ea;p(-ic i ; k (-q)i)r^;^(k',q')ak 1 (qi) a k 2 (q2) 

ea^(iatejqi)t')ea^(-iu^(q2)t')^^^ (60) 
The reduced form of the six-point function may be written as, 

I = A k (q)ezp(-KJ k (-q)t)r£;gj(k , q ) exp(i(J kl {c[x)t ) exp(- (ci 2 )t )A^(-c{ )exp(i^< (q')t' ) 

(a k (-q)a ki (qi)a k2 (q 2 )a k » (q" )) 

This may be simplified to 

/ = exp(—i w k (— q) t)exp(i w k (— q) £ )exp{—i w k " (q ) £ )exp(i w k " (q ) t )A k (q) 

V(-q")r k , k ^"(k',q) (61) 

Further, since the commutation rules, 

[ C k+q/2 C k-q/2 7 C k ' +q ' /2 C k '_ q '/ 2 ] = C k+q/2 C k' -q72^k-q/2,k'+q72 _ C k ' +q 7 2 C k-q/2 ( 5k+q/2,k' -q72 ( 62 ) 

This means, that the coefficents also have to satisfy, 

= V+?/2W + 9 )^k-q/2,k' +q /2^k 2 ,k' +q/2^q 2 ,-q-q - ^k' -q/2\1 + 9 )^k+q/2,k' -q' /2^k 2 ,k' -q/2^q 2 -q-q' (^3) 

First notice that A k (q) is either zero or one. Notice that Eq.( |5^), Eq.( ^) and Eq.( ^) are such that they 
suggest to us, 

r £; k2 ( k ,q) = unless A kl (- qi ) = 1 and A k2 (-q 2 ) = 1 (64) 

r k^k2( k ' q ) = (5 k 1 -q72,k 2 -q 2 /2^k 1 +q72,k+q/2'5k-q/2,k 2 +q 2 /2 
~ ^ kl +q 1 /2,k 2 +q 2 /2^ kl -q 1 /2,k-q/2^k+q/2,k 2 -q 2 /2 for A kl (-qi)A k2 (-q 2 ) = 1 (65) 

and 

r^;S(k,q) = for A kl (- qi )A k2 (-q 2 ) = (66) 

From this one may verify that the formula for the kinetic energy operator written down previously, Eq.( |57j ) is 
correct. Lastly, we must make sure that the following additional commutation rule holds. 

rK(k,q) r k - k *(k\q)-rK(k,q) I^£(k\q) 

= r£i;g(k' + q/2, q + q')5 k _ q/2 , k ' +q 7 2 - r£;£(k + q/2, q + q )<5 kW/a . k+q/2 (67) 

Verification of the above rule completes the confirmation that the formula for c k+q / 2 c k-q/2 m terms of the bose 
fields is exact. The algebra here goes through without a hitch expect in this respect, namely we have to assume 
that the A k (-q) = 1. 

To summarise let us write down the following: 
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THEOREM: 



t N N t 

c k+q/2 Ck -q/ 2 = n F (k)—5^ + (W T^y)[A k (q)ak(-q) + A k (-q)a k (q)] 

+ E E r k;;S(k,q)a kl (qi)«k 2 (q 2 ) (68) 
ki,k 2 qi,q2 

together with 

Tk^kaC 4 '^) = Aki(-qi)Ak 2 (-q2)[^ki-qi/2,k 2 -q 2 /2^ki+qi/2,k+q/2^k-q/2,k 2 +q2/2 

- ^k 1 +q 1 /2,k 2 +q 2 /2<5k 1 -q 1 /2,k-q/2'5k+q/2,k 2 -q 2 /2] (69) 

is an exact transformation of products of Fermi fields into Bose fields. Namely, all dynamical moments of 
c k:+q/2 c k-q/2 are recovered exactly, provided we identify the filled Fermi sea with the Bose vacuum. 

a k (q)\FS) = (70) 

and the above is also an operator identity, namely they satisfy proper mutual commutation rules. 

This may also be written in a simplified form as, 



t N N t 

c k+ q /2 Ck -q/ 2 = n F( k )j]^ 6 <i,o + (W T^y)[Ak(q)ak(-q) + A k (-q)a k (q)] 

+ E Ak+q/2- qi /2(-qi)A k _ qi/ 2(q - qi)a k+q /2- qi / 2 (qi) a k-qi/2(-q + qi) 
qi 

- E A k-q/2+ qi /2 (-qi ) A k+qi/2 (q - qi )a k - q / 2+qi /2 (<H K+qi /2 (~q + <Il) (71 ) 

qi 

The -Mr is needed to ensure that, X)k c iL c k = N (and not J2k c L c ^ = (AO) The (./ -Mr) is needed to ensure 



that the commutation amongst the c k+q / 2 c k- q /2 come out right. However, there are other finer points that need 
to be discussed. For example, do the thermodynamic expectation values also come out right ? That is, assuming 
that the temperature is finite, can we compute the dynamical density correlation function and show that it agrees 
with the free theory ? The finite temp, case is most unusual. It seems that for dealing with fermions at finite 
temperature, we have to assume that the Bose like-excitations are always at zero temperature with zero chemical 
potential, only the coefficients in the expansion acquire finite temp, values. All this is very strange and on shaky 
ground. Thus, I shall relegate the finite temp, case to the last section where the skeptical reader may have a field 
day demolishing my ideas. 



3.1 The Fermi Field Operator 

Also how about expressing the Fermi fields themselves in terms of the Fermi sea displacements ? For this one has 
to first construct a canonical conjugate of the density and then use it in the DPVA ansatz Eq.( ||) jlTJ] to compute 
the phase functional $. I have tried to construct a canonical conjugate of p ( any canonical conjugate suffices ) 
but it has not been successful. Just for future reference, it is useful to write down a formula for the total current 
operator. 

J - = E^^k = E E r k^ 2 « kl (qi)«k 2 (q 2 ) (72) 

k ki,k 2 qi,q 2 

= Ak 1 (-qi)A k2 (-q 2 )(^)5 kl , k2 J qi , qa (73) 



Therefore, 



(^)4ck - E Ak 1 (-q 1 )(^)at i (q 1 )a kl (qi) (74) 



m e — ' m 



k ki, q i 
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Thus an alternative but less systematic approach seems to be in order. For this choose, 

V>(x) = exp(-i U (x))exp(-i Ui(x))exp(-i f7 2 (x))(p(x))5 (75) 
here Ua(x) is a c- number and, 

E/i(x) =5Z[ak(q)/(k,q;x)+4(q)r(k,q;x)] (76) 

k,q 

U 2 (*) = E ^ r 2(ki,qi;k 2 ,q2;x)aj ci (q 1 )a k2 (q 2 ) (77) 

ki,qi k 2 ,q 2 

Now let us evolve this using the kinetic energy operator, 

V>(x,t) = e l * K V(x)e" 4 * K = exp(-i U (x))exp(-i f/i(x,i))eo;p(-i f7 2 (x,i))(p(x,i))2 (78) 

here, 

C/!(x,i) - ^K(q)e- 1 * ^>/(k,q;x) +a[(q)e 4 * ^/*(k,q;x)] (79) 

k,q 

E/ 2 (x,t) = ]T E ^(k 1 ,q 1 ;k 2 ,q 2 ;x)aj ii (q 1 )a k2 (q 2 ) e t 4 (^(n)-^^)) (80) 

ki,qi k 2 ,q2 

The basic philosophy it seems, involves observing that the commutation rules satisfied by number conserving 
objects such as c k+q/ / 2 c k-q/2 amongst themselves is independent of the underlying statistics. Thus one is lead to 
postulate convenient forms for such objects in terms of other Bosc fields. But not all approaches are likely to be 
correct. Take for example, the straightforward substitution, 

C !+q/2 C k-q/2 = fo I +q /2 6 k-q/2 (81) 

where c k are fermions and & k are bosons. This no doubt makes the commutation rules come out right but is quite 
obviously wrong because all the correlation functions are wrongly represented. 



4 Luttinger and Fermi Liquids 

Consider an interaction of the type, 



Hl = E ^Et A k(qK(-q) + A k (-q)4(q)][A k ,(-q)<v(q) + A k , (qjo*, (-q)] (82) 
q#o k , k ' 



. d 



gal(q) = c k (qK(q) + (^)A k (-q) £[A k , (-q)<, (q) + A k , (q)a£ (-q)] (83) 

k' 

The equations of motion for the Bose propagators read as, 

8 . u -t<r4(q)ot,(q)) 
(*^ - Wk(q)) ^) = <Vk"W^) 

, n v-",, , -*(Ta^(q)at,(q')) -^Ta^, (-q)^, (q')) 
+ (y)A k (-q) iJK" (-0) (rl) + A k " (q) <T1) ] (84) 

,.9 , -i<Ta k *(-q)4,(q')) 
(*-57+^(-q))- 



0t v (Tl) 

/ ^q v . , ,^ r . , -«(Ta[*,(-q)4(q')> . . -i(Ta k „(q)a k ,(q')) 
= -(^)A k (q)2^[A k «(q) ^ +A k «(-q) ] (85) 
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The boundary conditions on these propagators may be written down as, 

-^(-q^q)) - t (T4 (t -^(-q) a [,(q)) 
(Tl) (Tl) 

-»(raUq)<,(q')) _ ^{Ta^ (q)^, (q )) 
(Tl) (Tl) 



-i (3 uj. 
The boundary conditions imply that we may write 



-z(TaUq)<,(q')) ^ -i^q)^ (q')) 

= ^ esp(w„t) ^ (90) 



(Tl) 

= -(^)A k (q)2^[A k »(q) (Tl) + A k"(-q) (T1) ] (93) 



Define, 



-z(Ta^q)a^(q')) , 
2_, A k(-q) = Gi(q, k , q ; n) 



(86) 
(87) 



m = { 1 )E e^) (g9) 



i(Tat*(-q)at,(q')) -z(r4"(-q) t,(q')) 

= ^ exp(w„t) — — -s (91) 



(T1) ^ (T1) 

and, u„ = (2 7r n)//3. Thus, 

-»(TaC(q)qt,(q')) 6 k k ,^ , 
(«*. - - k (q)) ^ = — ]r 

^< V ^K < \ V^r a , ^ -»(ra^(q)a[,(q')) -»(Ta^ (-q) Q t, (q ')) 

+ (y)A k (-q) >JA k » (-q) k (yl) + A k » (q) — ] (92) 

^ , -i(Ta k "(-q)a k ,(q')) 
{iu n + Wk(-q))- 



(94) 



- ! (T4Vq)«I^ - . _ 

2^A k (q) (Tl) =G 2 (q,k,q;n) (95) 

Multiplying the above equations with A k (— q) and summing over k one arrives at simple formulas for Gi and G 2 . 



Gi(q,k ,q ;n) = A k -(-q) — 

i /3(iLJ n w k ' (q)) 

+ /„(q) [Gi (q, k' , q ; n) + G 2 (q, k' , q ; n)] (96) 

and, 

G 2 (q, k\ q ; n) = /* (-q) [d (q, k' , q ; n) + G 2 (q, k' , q ; n)] (97) 

and, 

/,(q) = (^)E f . Ak( " q i » (98) 
V ^ (iw„ - w k (q)) 
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In ID and 3D an analytical solution for /„(q) is possible. In ID we have, 

£ _ m i^n _|_ q ^ _ m itj n _|_ £ 

/n(«) = %(- — )(— )/n(— ^ — ) - 6(b - a) v„ ( — )(- — )ln( % -) (99) 

J»W Q\ q )\ 2 - K l V _ fc/ _ ELJ^n + 1> ^ > 9 V 27T A q> V a - + § ' V ' 

where, b = min(kf, kf — q) and a — max(—kf, —kf — q). The logarithm is to be interpreted as the principal value, 
that is, Im(ln(z)) — Arg(z) = tan~ 1 (y/x), x = Re(z),y — Im(z), and Re(ln(z)) — ln(r) and, r = (x 2 + y 2 ) 3 . 
The nice thing about this definition (not shared by non-principal value defns.) is that (ln(z))* — ln(z*). This 
is because — tt < Arg(z) = tan~ 1 (y/x) < ir. This may easily be solved and the final formulas for the Bose 
propagators may be written down as, 

G 2 (q,k',q';n) = G^q.k'.q'in) 

Gi(q,k',q';n) + G 2 (q,k',q';n) = Gi(q,k',q';n)/(1 - /*(-q)) 

Gi (q,k ,q ;n = — - — r ? ,. , „ 100 

-« /3 (1 - .f,t(-q) - /n(q))K - w k - (q)) 

n( ^' ' \ i 1 /n(-q)A k '(-q)3 q , q ' 

G 2 (q, k , q ; n ) = (—-5)7; 77? x f 1 YU - 7~v\ 101 

-%0 (1 - /*(-q) - / B (q))K - w k » (q)) 

r< c u' ' \^ /■ 1,' ' ^ c 1 ^ A k'(-q)^ q , q ' 
Gi(q, k , q ; n) + G 2 (q, k , q ; n) = (—7) 



i /3 7 (1 - .At(-q) - /n(q))(^„ - uv (q)) 
i<T<(q)a k ,(q')) «k,k'* q ,q' 



also, 



(Tl) -t /3(iw n - Wk(q)) 

+ l -z/3 n y j ^„-c k (q))(l-/*(-q)-/„(q))^„-^(q)) 1 J 

-z(T4"(-q)at , (q')> = 1 ^ Ah (q) A k ,(-q)g qiq , 



(Tl) /T v V '(iw„ + w k (-q)) (1 - /*(-q) - /„(q))(tw n - w k / (q)) 

The zero temperature correlation function of significance here is, 

-z<4,(q')a k (q)) (104) 

This may be obtained from the above formulas as, 

- Hh'W* - -(^.HtfVHO^ I ^7 (iM -„ t(q))(to - u ,.M )(1 - /;( - q) _;.(,)) « 

where G is the positively oriented contour that encloses the upper half plane. Thus the problem now boils down 
to computing all the zeros of (1 — f„(— q) — /«(q)) that have positive imaginary parts. In ID, we may proceed 
as follows, 

_ m iuj _|_ q 

1 - £("«) - /„(*) = 1 - <~){^-)H j _^* \ 



q 27T -fe/ 
+9(b-a)v q (-)(--)ln[ " 



2 



,m. , 1 . fc / 5 



g 2tt L -k f - nLJUL - 1 



1 wTO - 



2 

g 2 



+ f (6 " a K(^-)(-)/n[ _, J. t ] (106) 
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where a = max(—kf, —kf — q), b = min(kf, kf — q) and a = max(—kf, —kf + q), b — min(kf, kf + q). There 
are several regions of interest. 

(A) < q < kf This means, > —q > -kf or kf > kf - q > 0, —kf > -kf - q > -2kf, kf < kf + q < 2k f , 
— kf < —kf + q < Thus a = —kf, b = kf — q, a = —kf + q, b = kf. 9(b — a) = 1, 9(b — a ) = 1 Therefore, 



£ _ rn iui , q 

l-m-q)-fn(q) = l + v q (-)(-)ln[ " 2 



, 1 w m., ■ 
2n q 



q 2tt kf - 
-fc, _ HLJa _ 2 

■» q 2 - 

9 ^ 2 



1 TO (fcf +q/2) 2 + {^) 2 

i - «-,) - /.w = i + »,(^)( 7 )'-[ (t ;_ g/2)2+( j B)2 ] 



(^/ + g/2) 2 



'm_w *\2 
* 9 ' 1 



2 



9 

(kf + q/2) 2 + (^) 2 > T „ J 

u w ; v g ; = exv(-( q )( — )) 
(k f - q/2f + (^fY eXP[ ( to K V J 

C^ f = -[(kf + q/2) 2 (kf ,/2) 2 e, P (-(^)(l))]/[l - exp(-(^)(± 

y fit U fi lit U t 



-(^)(-) 
to w 9 



)(-))] 
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We want to find a root of this that has a positive imaginary part. 



w = i ( — ) 



\ 



(kf + - (fc/ - q/2) 2 exp(-e-^)(±)) 

l-exp(-( 2 -^L)(±)) 



(107) 

(108) 

(109) 

(110) 
(111) 

(112) 



The next case is, 

(B) -k f < q < 0. In this case, < k f + q < kf, -2 kf < -kf + q < -kf, k f > -q > 0, 2 kf > kf - q > kf, 
> —kf — q > —kf, a = —kf — q, b = kf, a = —kf, b = kf + q. Thus, 8(b — a) = 1, 6(b — a ) = 1 For case (B) 

,m w 1 N , ,~ kf 



9 g v 2tt' L -k f 



l-r n (-q)-f n (q) = l + v q (-)( — )ln 

q 



m ilu q 

9 2 - 

m iuj I q_ - 

9 ^ 2 



1 . ,771, 



+ W 9(^)(-)^ 



9 ^ 2 ■ 



9V 2vr" 9 '"'~ L fe/ - 
Thus the pole for both cases (A) and (B) is given by, 



m iu 9 - 

9 2 



(113) 



TO 



\ 



(fc/ + g/2) 2 - (fc/ - g /2)2 ea ;p(-(^)(^)) 

1 



(114) 



The quantity under the square root is positive in both cases. Third, we have the case, 

(C) kf < q < 2 kf , this means, —kf > —q > —2 kf, also, 2 fc/ < fe/ + q < 3 kf, < — fc/ + q < kf, 
> kf — q > — fc/, —2 fey > — kf — q > —3 fc/. a = — fc/, b = kf — q, a = —kf + q, b = kf. 9(b — a) = 1 and 
0(6 — a ) = 1 This leads to the same root Eq.( 114). In case 

(D) -2 kf < q < -kf, -kf < kf + q < 0, -3 kf < -k f + q < -2 kf , 2 kf > -q > kf , 3 fe/ > fe/ - q > 2 kf , 
fe/ > —kf — q > 0. a = —kf — q, b — kf , a = —kf, b — kf + q. 9(b— a) = 1 and 8(b — a ) = 1. In this case also we 
find the root given by Eq.( 114). (E) 2 fe/ < q < oo, —2 kf > —q > —oo, 3 fe/ < kf + q < oo, kf < —kf + q < oo, 
— kf > kf — q > — oo, —3 kf > — fe/ — q > —oo. a = —kf, b = kf — q, a — —kf + q b = kf. 8(b — a) = 0, 
0(6 — a ) = 0. This case is somewhat different. 



1 - /*(-<?) - /«(«) = 1 + v q (-)(^)ln[^- 

q Z7r —kf 



9 ^ 2 



m la; I q ■ 
9 + 2 



1G 



g zir 



(115) 



This when solved gives the same result as Eq.( 114). So too does the final case (F) which is 
(F) — oo < q < — 2 kf, — oo < kf + q < —kf, — oo < —kf + q < —3 kf, oo > —q>2 kf, co > kf — q > 3 kf, 
oo > —kf — q> kf. This case also leads to the same result namely, Eq.( 114). Therefore the final result may be 
written as, 



, t , \ ( )) = r I 1 Afc(-g)A fc '(-g)^^ 

^ k ' lq >a ^ q>> ^W<?) +u k {q)){u R {q) + ^ ( q ))(^)(^)2(^ R (q)(cosh(X( q )) - 1) 



Here we may write, 



m v a 



WR{q ) = (M) /( fc /+g/ 2 ) 2 - - qf2) 2 exp(-X(q)) 



1 - ezp(-A(g)) 



(116) 

(117) 
(118) 
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In other words, 



(4 c fe> = n F (k) + {2Trk f ) 



+00 



dqi Afe_ gi / 2 (— ?l) 



2tt 2w fl (gi)(a; fl (gi) + LJ k _ qi/2 ( qi )) 2 (^-)(cosh(X( qi )) - 1) 



_ (Jhrfc ) [ +co d Ji W^) (119) 

In order to see how good the present theory is, it is desirable to compare these results with the Calogero-Sutherland 
model or more specifically with the spin-spin correlation function of the Haldane-Shastry model. This is given by 
@, p = 1/2, a = 2, m = 1, k f = tt/2. 

(0\&(x,0)9(x ,0)|0) = / — exp(i k (x - x))(c{c k ) 
J-00 2vr 



^(r(i/2))2(r(i)) 2 



-1 



7T > 7T / 

exp(i — {x — x ) vi)exp(i — (x — x) V2) (120) 



(r(3/2)) 2 ^"(z-z')*™ 



(r(l/2))2(r(l)) 2 ^ ' (2n)! (2n+l)= 



For this we have to use the interaction given by (V(x) 



a(a — 1) \ 



v 9 = -a(a-l)(n\q\) (121) 



4.1 Qualitative Conclusions 

It is clear that one of the features of the Luttinger liquid that is a result of the Mattis-Leib solution is absent, 
namely the discontinuous dependence of the momentum distribution on the coupling strength as the latter goes 
to zero. This may therefore be an artifact of the Luttinger model and not a generic feature of all ID systems. Also 
the highly nonperturbative manner in which the momentum distribution of the interacting system approaches 
the noninteracting one as the coupling goes to zero may be seen quite easily. The other point that may be seen is 
that (pq^o) = indicating that there is no Wigner crystallization at any density. Furthermore, the vanishing of 
the quasiparticle residue is seen only for sufficiently large values of the repulsion between the fermions. Thus the 
conventional view that the fermi surface is destroyed for arbitrarily weak repulsion in case of Id fermions does 
not seem to hold up in this case. 



4.2 Momentum Distribution in 3D 

For this one has to compute the expressions in 3D: 



he ■ & I Q I 

2 u m m ,*Wn-Cq- ^ 



/n(q) = / dk k 



(2n)*J Q k\q\' ^ n _ £q + M 



777 tUJ 

-(-IL)6(k F -\q\-k)ln' 



q 

k\q\ 



£fc f + £k 



+ (rrrMkl + k - k M k f - ( k - W)H~ — ( 122 ) 
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This may be split into several cases 
(A) < |q| < k f 

/n(q) = 



(B) k f < \q\ < 2 k f 



/n(q) 



fe|q| 



S-(^){- / ' dk k Inip^ 1 



W\q\ Jl Jo 

kf-\q\ 



+ 



J 

JO 



lUJ n - e q 
fe|q| 



fe|q| 



dk k Jn(-= l—fir) 



fc F -|q| 8Wn-£q+ V 

<! / dkklnC n ~^ 

Jo 

dk k ln{ 



WW 

rkf 



Mil 



2W„ - £q + 



k\q\ 



M-kf 



lOJ n - e q + 



klg 



■)} 



(C) |g| > 2 fc, 



/n(q) = -7T^2(fr) / dkkln{- 



k\q\ 



^ n - £q + M 



(A) 



.TO 2 , , 1 . . 



/■ \2j A^n £q 

+ (iw„ - e q ) ln( 

+(^)(^n-e q )} 

/• \, /• n • / fcfM ■ \, / ■ x / kf\q\ . 

+ to (iuj n )ln(iLu n ) - m iu n - m(e„ h tuj n )ln(e„ h iLO n ) + m(e q h iw n )\ 

m m to 



(B) 



/n(q) 



m 



to 



.TO 2 , 1 . 

+ ( 2^ ){ 2 ( ^'' 



+(iu>„ - e q ) Zn(- 



,-W) 2 -i^ l + £q -W) 2 

TO 2 TO 



iw n + e q 



kf\q\ 



+ (W)(i w „-e q )} 

TO 

m /. n • / kf\q\ . ,, , fcflgl . . . fcflgl 
+ to (iw„)in(iw„) - to iu n - TO(e q h iw„)/n(e q h iu n ) + TO(e q — h iu) n )\ 



TO. 



TO. 



TO 
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(C) 



iio n - e q + 



fc/kl 



fc/kl 



m' 
2k" 



2 {g^n 



, fc/kl 



TO 



+ (iu n - e q ) 2 Zn(- 



) z — (zw„ -e q )}] 



zw„ - e a - 



fc/kl 



TO 



(128) 



5 Fermions at Finite Temperature 

It should be clear from the previous sections that fermions at finite temp, is a major headache. A straightforward 
generalisation is not working out for unknown reasons. In this section we shall not distinguish between N and 
(N), things are complicated enough as it is ! Take for example, the expectation value of the number density at 
finite temerature, 

(4ck) = n F (k) + A k-q 1 /2(-qi)(aj c -q 1 /2( c li) a k-q 1 /2(qi)> 
qi 

- A k+ qi /2(-qi)(aI +qi/2 (qiK+ qi /2(qi)} (129) 

qi 

Assume the n F (k) are all evaluated at zero temp, as in the previous sections. Now, taken at face value, one 
is obliged to compute the thermodynamic expectation values of the Bose occupation probabilities assuming the 
chemical potential for the Bosons is zero (this means that we are allowed to create and destroy any number of 
Bosons) . Now such a calculation yields an infinite answer for (c k c k ) as the sum over all qi diverges (is proportional 
to the total number of fermions). This is the reason why a more nefarious approach may be necessary in dealing 
with fermions at finite temperature. For this the only guide is that all the finite temp, dynamical correlation 
functions involving the number conserving object c k+q/ / 2 c k-q/2 should be correctly reproduced. Hopefully the 
commutation rules involving these number conserving products are not damaged in the bargain. For this let us 
start with the simplest case namely, (c k+q//2 c k _ q / 2 ) ■ We know what the answer should be, that is, 

(4+ q /2 c k-q/2> = Sq=on F ^(k) (130) 

where, 

n Ft0 {k) = [exp(/3(e k - fx)) + l]" 1 (131) 
For this to happen we have to do the following, 

(1) Treat the bosons as before assuming that they are always at zero temperature and with zero chemical potential. 

(2) Assume that all finite temperature effects are lumped into the coefficients. 

(3) Fix the coefficients as before by demanding that the finite temp, dynamical moments of the number-conserving 
products come out right. 

For the dynamical four-point function to come out right we must ensure that, 

A£(q) - yW^k + q/2)(l - n Ffj (k - q/2)) (132) 

Unfortunately, the coefficient A^(q) is no longer either zero or one. This renders matters even more difficult. But 
the dynamical six-point function is still the same and has to be correctly given as in the zero-temp case. 

1 = ( C k+q/2 C k-q/2C^ +q / /2 C k '_ q72 cj c „ +q „ /2 C k "_ q » /2 ) 

= [(1 - n F ,/3(k - q/2))(l - n Ffj {k - q'/2))n F , /3 (k + q/2)^ k+q/2ik "_ q » /2 ^ k _ q/2jk ' +q ' /2 (5 k '_ q ' /2ik » +q " /2 
- (1 - n Ffj (k - q/2))n Ffj (k + q /2)n FJj (k + q/2)^ k _ q/2;k » +q " /2 (5 k+q/2;k /_ q72 ^ k / +q72 k »_ q " /2 ] (133) 
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In terms of the Bose fields we have, 



/ = A£(q)Af,,(-q")r w q w ,? (k',q) (134) 



Equating these two we arrive at a formula for the coefficient T. It is somewhat different from the zero temp case. 
But we can remedy that by choosing, 



r££(k,q) = F(A^ i (-q 1 )A^(-q 2 ))[ v /(l - + qi /2))(l - n Fj/3 (k 2 + q 2 /2)) 

<5ki-qi/2,k 2 -q2/2e>ki+qi/2,k+q/2<!>k-q/2,k2+q2/2 



'n Fi(3 (ki - qi/2)n F ,/3(k 2 - q 2 /2) 

<5 k 1 +q 1 /2,k 2 +q2/2 (5 k 1 -q 1 /2,k-q/2^k+q/2,k2-q2/2] (135) 

and, 

F(x) = 1 if x ^ and F{x) = for x = (136) 

For any temperature above zero F(x) = 1 always. But for exactly at zero temp, the original case is recovered. 
But more importantly, all the right sort of dynamical correlation functions are recovered if we choose in addition, 
the following formula for the kinetic energy operator, 

X = ]>> k (q)4(q)a k (q) (137) 

k,q 



and, 

^(q) = F(Af(-q))(^) (1™) 



m 



These choices ensure that all the finite temp, dynamical moments of c k+q / 2 c k-q/2 are correctly recovered. 

5.1 Demonstration of Emergence of Superconductivity 

Consider two interactions, one the coulomb repulsion and the other electron-phonon interaction : 

HI = 52W E[ A k(<lH(-q) + A£(-q)4(q)][A£,(-qK,(q) + A* (qjo*, (-q)] (139) 
q#o k,k' 

H e - phon = £ ^ Et A k(q)«k(-q) + A£(-q)4(q)][6 q + &L q ] (140) 

q#0 V 2 k 

The displacement operators are, 

= & q + &-q ( 141 ) 

and, 

* q = (!)(&_<, -&t) (142) 
The equations of motion for the Bose fields are, 

i |a* k (q) = ^ k (q)a k (q) + (^)Af(-q) £[A£ (-q)a£, (q) + A* (qjaj (-q)] 

k' 

+ ^A£(-q)[& q + &! q ] (143) 
= TrElV ^K' ("<0 + A k ' (-q)a^ (q)] + ^P_ q (144) 

k 



2 n c 

dt 



f «x p q = — "*-q ( 145 ) 
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* ^ a k(-q) = -"k(-q)o£(-q) - (f )A^(q)^[Af,(-q)4,(q) + A* (q)a£ (-q)] 

k' 

-^(q)[6 q + 6L q ] 

1/ 2 



- Wk(q)) = 5 k , k 'tf q , q '5(t) 

« a „ „ „ -i(T a?„(q)oL(q')) fl 4*"(-q)4'(q')) 
+(^)a£(-q) £[a£„ (-q) + <" (q) /n ] 



V kV ^f^ L k v M; (Tl) k (Tl) 

+ ^r ki qj (ti) 



(«^+«k(-q)) ^ = 

-(ir)A^(q) 2JA£„ (-q) k , , k + A£„(q) * k - 



M "A^ n r i<r ^^ ' (q)> 



The corresponding propagators with phonons are, 



. a -%(T Pq 4,(q')) _ 20, -i<T Ai q o+,(q')) 



at (Tl) v i ' (Tl) 

,a -«{T^ q flt,(q)) _ j tf _ q ^ ff -i(r^(q)at( q ')) -»(T a^(-q)a 

dt (ti) ~ l J2J A k"i ^ (ti) +VW (T1) 



k 

+( i £ a ) ' ' ' 1 k 



O q ,-^P* at (q')> 



2 ' (Tl) 

Transform to frequency domain and we have the corresponding equations : 

Gi(q,k ,q ;n) = A k <(-q) — 



-i /3(iw n - w k '(q)) 

+ /n (q) [Gi (q, k' , q' ; n) + G 2 (q, k' , q' ; n)] + g n (q)G 3 (q, k' , q' ; n) 

G 2 (q,k',q';n) = /*(-q)[Gi(q, k', q'; n) + G 2 (q, k' , q'; n)\ + .g*(-q)G 3 (q, k' , q' ; n) 

, , 2 / / 

i u n G 3 (q,k ,q ;n) = ( — — )G 4 (q,k ,q ;n) 

i w„ G 4 (q,k ,q ;n) = i [Gi(q,k ,q ;n) + G 2 (q,k ,q ;n)] + (-^)G 3 (q,k ,q ;n) 

The solutions may be written down as, 

/ / 2 ii ii 

G 3 (q,k ,q ;n) = -( q 2 )[Gi(q,k ,q ;n) + G 2 (q,k ,q ;n)] 
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G 2 (q,k',q';n) - ( ^" (q) )Gi(q,k',q ;n) (156) 

i J- n iqj 



rn(q) = /:(-q)- ? |^< ? :(-q) (157) 



and finally, 



and, 



G 1 (q, k' , q ; „) = <5 q q , F ^ )K ' { q) (158) 
qq -t 0(i uj n - w k /(q)) 



here. 



PCn^-M /»(q) , 2 q g n (q) x 

^ l(q) - [1 _ T^(q) + ( ?fT^ )( T^ryq) )] 

(Af(-q)) 



(159) 



V ^ {iui n - wk(qj) 



k 



, x ,|M q | 2 ^ (A^-q)) 2 MRn 
1/ ^ (iw„ - Wk(qj) 

For the four-point function we have to proceed as follows, the quantity of interest is(with phonons and coulomb), 



2(n«L 1 1 (qi)«L 2 2 (q 2 )a^( q ;)a^( q ' 1 ))) 



for this we as usual decompose as follows, 



4 (T« (qi) < (q 2 )a]> (q 2 )aj> (qi))) 



= ^ exp( Wn t 1 )i(T(aC 1 (qi) <(q 2 )a^( q ;)a^( q ;))) (162) 

n 

(i uj n - cj kl ( qi ))^<T« ( qi )a& (qa)a^ (q 2 )ajj (q' x ))) 
= Fx(kx q! n,k 2 q 2 i 2 ,k 2 q 2 t 2 ,ki q' x 4) 

x [Fi(qin,k 2 q 2 i 2 ,k 2 q 2 t 2 ,ki qi t[) + F 2 (qm, k 2 q 2 i 2 ,k 2 q 2 i 2 ,ki qi ti)]/(l - r n (qi) - T^-qi)) (163) 



fi(qin,k 2 q 2 t 2 ,k 2 q 2 tj,^ qj i x ) 

E ^(-q 1 ) F 1 (k iqi n,k 2 q2 i 2 ,k; q ; 4ki q ; 4) (164) 

^ (iw„ - wki(qij) 



ki 



F 2 (q!n,k 2 q 2 i 2 ,k 2 q 2 t 2 ,^ qj ^) 

^(qi) 
(ioj n + w kl (-qi))" 



X! r — yr^2(kiqin,k 2 q 2 t 2 ,k 2 q 2 t 2 ,k'i qi t[) (165) 
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5.2 Bosons at zero temperature 

In this subsection we do what we did earlier except that here we are bosonizing the bosons. Let 6k and 6^ be the 
Bose fields in question. Analogous to Fermi sea displacements we introduce Bose condensate displacements. 

6 k+q/2 fe k-q/2 = N <5 k , <Sq,0 + (VN) [5 k +q/2,0 dk(-<j) + ^k-q/2,0 4 fa)] 
+ Tl ( k ' q ' t ll)4+q/2-q 1 /2( t ll) d k-q 1 /2(-q + Ql) 

qi 

q/2+q 1 /2( C ll) d k+q 1 /2(-q+ Ql) (166) 

qi 

To make all the dynamical moments of 6j c+q ^ 2 6k-q/2 come out right and the commutation rules amongst them 
also come out right provided we choose 2\ and T 2 such that, 

(V7V) 2 5 k+q/2 , V'-q'72,0^ Tl ( k '' q '' qi )^ k ^ q )4'W/2-qi/2( qi ) d k'- qi /2(-q' + <U.)<4" fa")) 

qi 

+(^) 2 W2,oVw72,oE T2 ( k ^ q '' q i)< dk ^^ 

qi 

N < 5k+q/2,0'5k+q/2,k"-q"/2^k-q/2,k'+q72^k'-q72,k"+q"/2 (167) 

This means, 

Ti(-q-q72,q',-q) = 1; q £ 0; qV (168) 

r 2 (-q72,q,-q)=0; q^O; qV0 (169) 
in order for the kinetic energy operator to have the form, 

* = $>k d{ 1/2)k (k) d ( i /a)k (k) (170) 

k 

we must have, 

£k+ qi /2ri(k + qi/2,0,q 1 ) + e k _ qi/2 T 2 (k- qi /2,0, qi ) = 5 k , qi/2 e qi (171) 

In order for, 

k 

we must have, 

Ti(k + qi/2,0,qi)+T 2 (k-qi/2,0,qi) = (172) 

In order for the commutation rules amongst the 6 k+q / 2 6 k - q /2 to come out right, we must have in addition to all 
the above relations, 

T 2 (q/2,q,-q)=0 (173) 
T 2 (-q/2,q,q + q')=0 (174) 



Any choice of T\ and T 2 consistent with the above relations should suffice. From the above relations (Eq.( |171| ) 
and Eq.( |l72|)), we have quite unambiguously, 



T 1 (k + q 1 /2,0,q 1 )=5 k)qi/2 (175) 

T 2 (k- qi /2,0, qi ) = -5 k , qi/2 (176) 

T 1 (-q-q72,q,-q) = l (177) 

and, 

Ti(k, q, qi) = 0; otherwise (178) 

T 2 (k, q, qi) = 0; otherwise (179) 
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This means that we may rewrite the formula for &j c+q/ / 2 ^k- q /2 as follows, 

6 k+q/2 5 k-q/2 = MSk,0<*q,0 + (V^V) [<5k+q/2,(A(-q) + ^k-q/2,o4( ( l) 

+ d (i/2)k+ q /2( k + q/ 2 ) d (i/2)k-q/ 2 (k - q/2) 

E< l/2 (qi) d qi /2(qi) ( 18 °) 

qi 

Consider an interaction of the type, 



or, 



H I = (yiyj^yit^+qAO dk(-q) +5 k _ q / 2 ,o <4(<j)Pk'-q/2,0 d k'(q) + ^k'+q/2,0 

4(-q)] (181) 

q#o k ,k' 

^ = (f)£^-q/2(-q) +< /2 (q)]K/ 2 (q) +d t _ q/2 (-q)] (182) 

q#0 

and the free case is given by, 

^0 = £ £ q rf (l/2)q( C l) d (l/2)q( f l) ( 183 ) 

q 

The full hamiltonian may be diagonaliscd as follows, 



H = 5>q4A (184) 



and, 



/q = ( W q/2 (q) + ( ^-7T^ t -) 2d -q/2(-<l) (I 85 ) 

Z CJq Z CJq 

/lq - ( -^ + ^ + ^ )^ q/2(q ) + (^±Ma)^t ( _ q) (186) 

^q/ 2 (q) - r q+ ! q + ^ q )^q - ( -^ + ^ + ^ )*/t q (187) 
Z CJq Z CJq 

,t_ q/2( _ q) = ^q + eq + ^q ^^ _ ( ^q + Cq + ^ } 1 ^ ^ 

CJ q = y / e q + 2 / 9 Wq£q (189) 

From this we may deduce, 

(< 1/2 )q(q)rf(l/2)q(q)) = Z^+^±^q (190) 

From this it is possible to write down the filling fraction. 
FILLING FRACTION 

f = N /N = 1 - (l/N) E( rf (l/2)q(q)^(l/2)q(q)) (191) 

q 

or, ^ 

/o = AViv = 1 - (l/ 27 r 2 po) r dq f ^ + ^ + P^ (ig2) 

Here the lower limit is necessary to ensure that the state 

rf(i/2) q (q)|G) 

has positive norm, and the value of q m in is given by, 

_ + e <? m ,„ + P0Vq min = (193) 
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First define, 

f UJq + 6 q + PQVq^l 

2 u, 



Aq = fq ■ .q , ^ q )j (1!)4) 

g q = ( -^ + 2 ^ + P0 ^ )* (195) 

The density operator is, 

p q (t) = >/jV[d_ (1/2)q (-q)(t) + d|i/2)q(<0(*)] + E d (i/2)k+ q /2( k + q/2)(t)d ( i/2)k-q/a(k - q/2)(t) (196) 

k 

d_ ( i /2)q (-q)(t) - A q .f_ qe -* - Bq/te* **' (197) 

4l/2)q(q)W = ^ - Sq/-qe- < (198) 

d(i/2)k-q/2(k - q/2)(t) = A k _ q/2 / k _ q/2 e- 4 ""■-«/»« - Bk- q / 2 /I k+q/2 e 4 Wk -<* /2 ' (199) 
d (i/2)k +q /2( k + ^ 2 )(*) = A k + q/ 2 / 1 [ +q/2 e l - Bk +q /2/-k- q /2e-' ^/ 2t (200) 

Now define, 

S>(qt) - (p q (t)p_ q (0)) - iV<[d_ (1/2)q (-q)(t) + d| 1/2)q (q)(t)][d (1/2)q (q)(0) + d^^-qXO)]) 

E( d (i/2)k +q / 2 ( k + q/2)(t)d (1/2)k _ q/2 (k - q/2)(i)dj 1/2)k ,_ q/2 (k' - q/2)(0)d (1/2)k / +q/2 (k' + q/2)(0)) 
k,k' 

= N(— ) ea;p(-i w q i) 

+ E(^ +q /2/- k -V2^' Wk+q/2 '4- q / 2 / k - q /2e-'" k -^ i 4^ q/ J^ q/ /, +q/2 /!,_ q/2 ) 

k,k' 

+ E( B H q /2/- k -V2^ Mk+ ^ i S k - q /2/! k+q/ /^-' 2t 5 k '_ q/2 /_ k ^ q/2 B k ' +q/2 /l^ 
k,k' 

S>(q,t) = iV(— ) exp(-i w q i) 

Wq 

+ ^ea:p(-j (w k+q/2 + w k _ q/2 )t)[B k+q/2 A k _ q/2 + B k+q/2 B_ k+q/2J 4 k _ q/2 A k+q/2 ] 

k 

5 ,< (q,t) = A^(— ) exp(i w q i) + ^ ea;p(i (w k+q/2 + ^k- q / 2 )i)[S k _ q/2 A k+q/2 + B k _ q / 2 B_ k _ q/2 A k+q/2 ,4 k _ q/2 ] 



and, 

5>(q,t) = (p q (t)p_ q (0)) 
S<(q,t) = <p_ q (0)p q (t)) 

From the above equations, it is easy to see that there is a coherent part corresponding to the Bogoliubov spectrum 
and an incoherent part which is due to correlations and is responsible (hopefully) for the roton minimum. 
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